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Abstract 

ly-j ■ We investigate ground state configurations for a general finite number N 

of particles of the Heitmann-Radin sticky disc pair potential model in two 
p | | dimensions. Exact energy minimizers are shown to exhibit large microscopic 

fluctuations about the asymptotic Wulff shape which is a regular hexagon: 
There are arbitrarily large N with ground state configurations deviating from 
the nearest regular hexagon by a number of ~ iV 3 / 4 particles. We also prove 
that for any N and any ground state configuration this deviation is bounded 
above by ~ iV 3 / 4 . As a consequence we obtain an exact scaling law for the 
fluctuations about the asymptotic Wulff shape. In particular, our results give 
a sharp rate of convergence to the limiting Wulff shape. 
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Jg 1 Introduction 

q \ It is a fundamental problem in statistical and solid mechanics to derive material 
properties from atomistic interaction models. One of the important open questions 
is to explain theoretically why atoms at low temperature arrange on periodic lattices. 
In fact, at very low temperature, a system of atoms is not only observed to arrange on 
the atomic length-scale. The atoms are moreover seen to form large clusters whose 
overall polyhedral shape is given by the Wulff shape which is obtained through a 
surface energy minimization problem of an effective continuum model. 

The first crystallization result was obtained by Heitmann and Radin in [3]: For 
a very special pair interaction model in two dimensions, see ([T]) and §Z§ below, 
they showed that ground state configurations at zero temperature are subsets of a 
triangular lattice with nearest-neighbour distance 1. (See also [3] for an equivalent 
problem in a different context.) The overall shape of energy minimizers for this 
model was analyzed in [TJ: We proved that in the large N limit, the unique rescaled 
asymptotic shape of ground states is given by a regular hexagon, which forms the 
Wulff shape of a limiting continuum variational problem. 
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In this paper we address the problem of characterizing the discrete energy min- 
imizers for finite N. In particular, we determine an exact scaling law for the mi- 
croscopic fluctuations about the limiting Wulff shape: Ground state configurations 
deviate from the nearest regular hexagonal shape by at most ~ iV 3 / 4 particles. This 
bound is sharp as there are highly degenerate discrete ground states for infinitely 
many N, which deviate form the nearest regular hexagon by ~ iV 3 / 4 particles. It 
thus turns out that there are surprisingly large fluctuations even at zero temper- 
ature, as rearrangements of surface atoms would only lead to deviations of order 
~ N 1 ' 2 . 

Our model energy, the Heitmann-Radin sticky disc model, assigns to N particles 
in two dimensions with positions x±, . . . ,x^ EM. 2 the pair potential energy 

E(xi, ...,x N ) = y]VgB.(\xi - Xj\), (1) 
with interaction potential given by 

{+00, < r < 1, 
-1, r = l, (2) 
0, r > 1. 

Such a potential can by interpreted as a very rough approximation of a hard core 
Lennard- Jones-type potential with infinitely short interaction length, see Figure [D 
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Figure 1: The Heitmann-Radin sticky disc potential. 

Based on the crystallization result in |3J H], for this model it was shown pQ that 
ground states approach a unique hexagonal shape, in the sense that the re-scaled 
empirical measure of any iV-particle minimizer [x\ , . . . , 

N 

[l N = iV" 1 ^ S N -1/2 X (N) , (3) 
3=1 

satisfies, up to translation and rotation, 

A*jv A poXh, (4) 
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where Xh is the characteristic function of the regular hexagon h with corners ±^= , 

273(^3)' = ' = 2v^(v^)' an< ^ ^° = 73 * s ^ e density °f atoms in a triangular lattice 
with nearest-neighbour distance 1. 

For finite N we will measure deviations \\fiN — PoXh\\ from the nearest regular 
hexagon h in terms of a suitable norm, see ([7j). Our main results are as follows. 

• A simple observation shows that there always exists a ground state configura- 
tion which agrees with a suitably chosen regular hexagon up to a surface con- 
tribution of ~ N 1 / 2 atoms, i.e., \\[xn — PoXh\\ < iV -1 / 2 , see Figures Eta),(b),(c) 
and El 

• For many arbitrarily large N there exist ground states which are highly de- 
generate: they differ from the nearest regular hexagon by ~ iV 3 / 4 atoms, i.e., 
\\pN — PoXhW > ciV _1//4 , see Figure [2]^d) (and Theorem 12. 6p . 

• Finally, this rate is sharp, in the sense that \\{an — PoXh\\ — CN~ 1 / 4 for all N 
and any ground state (see Theorem 12.41) . 

It is interesting to note that there are also infinitely many N for which the ground 
state is (up to rotation and translation) unique and for which the ground states of 
iV + 1 particles are highly degenerate. 
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(a) (b) (c) (d) 

Figure 2: (a) The unique (up to rotations and translations) ground state for 61 
atoms, cf . Proposition 12.51 (b) A ground state for 62 atoms, cf . Proposition 12.71 (c) 
Another ground state for 62 atoms obtained from (b) by rearrangement of surface 
atoms. New positions are indicated by white lattice sites, (d) Yet another ground 
state for 62 atoms obtained from (c) by rearrangement of a large surface layer of 
atoms. 

As an offspring, we obtain an alternative, perhaps conceptually less transparent 
but more elementary, proof of our result in [1] of convergence to the hexagonal Wulff 
shape, which also provides a sharp rate of convergence. 

We also remark that obviously our results hold true for any pair interaction 
model with interaction range less than the next nearest neighbor distance y/3, for 
which ground states are known to crystallize on the triangular lattice. This is e.g. 
the case for the Radin soft disc potential Vr, for which V&(r) = 24r — 25 on [1, ||] 
and which equals Vhr elsewhere, see [5], and thus for any V with Vr < V < Vrr- 
However, it seems that there are no examples known which allow for an elastic 
range, i.e., V quadratic in a neighborhood of the equilibrium point 1, as e.g. the 
short-range soft potentials considered in jTJ Section 2], although we conjecture our 
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results also to remain valid for these models. In fact, we expect that even long-range 
Lennard- Jones type interactions lead to Wulff shapes in the infinite particle limit. 
Theil remarkably succeded in proving crystallization under crystalline or periodic 
boundary conditions for a class of such models [7]. In our context, however, such a 
result would have to be proven, at least asymptotically, for the free boundary value 
problem, which appears to be a very difficult problem. 

Our results explain, in a model problem, fluctuations of atomic configurations 
at zero temperature. The rapid change of unique and highly degenerate ground 
states as N varies leads to an interesting discrepancy between the canonical and 
grandcanonical Gibbs ensemble restricted to absolute minimizers of the potential 
energy. It would be interesting to explore to what extend these results persist at 
low but finite temperature. In this regime the canonical Gibbs ensemble is expected 
to behave more regularly since even for particle numbers N with unique ground 
states a construction as in the proof of Theorem 12.61 shows that the second lowest 
energy level is highly degenerate. In line with our analysis at zero temperature, we 
expect to observe surprisingly large entropic effects for the statistical fluctuations 
about the hexagonal Wulff also at finite temperature. 

2 Analysis of ground states 

Let Sn = {xi, . . . ,xn} denote a general configuration of atomic positions Xi G M 2 
with finite energy. Hence \xi — Xj\ > 1 for all % ^ j. By viewing the atomic positions 
as vertices that are linked to each other by an atomic bond whenever their distance 
is equal to 1, we have assigned a graph to every configuration of atoms, whose set 
of vertices is Sn and whose set of edges B consists of their atomic bonds. 

An element of B is called a boundary edge and said to belong to dB if there are 
less than two vertices that are connected with both of its endpoints. We call those 
x G Sn for which #{y G Sn '■ \y — x| = l}<5 boundary points and denote the set 
of boundary points by 8Sn- 

The main result of Heitmann and Radin in [I] states that ground states are 
subsets of a triangular lattice with minimal interatomic distance 1. Of course any 
rigid motion maps ground states into ground states, and we will therefore without 
loss of generality only consider atomic configurations Sn which are subsets of the 
lattice 



For later reference we include a precise statement of the Heitmann- Radin results. 

Theorem 2.1 ([2111]). Let Sn be a minimal energy configuration of N > 3 atoms 
with set of atomic bonds B. 



(ii) Up to a rigid displacement the elements of dB form a simple closed polygon P 
with vertices on C Sn consists precisely of all the points inside and on P. 





(i) Then = [3N - y/12N - 3 J . 



(in) #dS N = #dB 



[3 - - 3J . 
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As an immediate consequence of (i) we obtain that the ground state energy for 
N atoms is given by 



-2#B = -2 



3N - V12N - 3 = -6iV + 2 Vl2N - 3 



(6) 



We are interested in the ground states for large finite N. So in the following 
results it is no loss of generality to (tacitly) assume that N be sufficiently large. 

Let Sn be a ground state. We choose an ordering of the elements of dSw = Sn^P 
by running through P counter-clockwise as v±, v 2 , ■ ■ ■ , v m , v m+ i = v 1 . For v G dS^ 
we denote by <p v the interior angle of P at v (taking values in {|, ir, y }) and set 

fi = <P Vi =^~{v i+ i - Vi, - Vi). 

Lemma 2.2. Let Sn be a ground state. 

(i) Then tfi = -^-,<fi + i = . . . = ipi + j = n implies that <^i+j+i 7^ -y- (Here and in 
the sequel all indices are understood modulo m.) 

(11) Iftpi = f, then <pj y£ f for j £ 1}. 

(Hi) There is at most one i G {1, . . . , m} such that ipi = |. 

(iv) If fi = |, then <pi-i,<p i+ i G {n, and not both of these angles are equal to 
n. 



Proof. Note first that for any configuration of atoms filling the lattice inside and 
on some polygon P whose segments are aligned with the nearest neighbor bonds 
in £, there always exists a boundary vertex v of degree deg(t> ) < 3 which lies in a 
prescribed half plane that intersects P. (Simply take a suitable corner of P.) 

(i) If this were not the case, i.e., we had a configuration as sketched in Figure El 

Vi+j+2 Vi-x 
• oooooo» 



V i+j + l 



Figure 3: Free lattice slots between atom and v i+ j +2 . 

we could successively move atoms of degree at most 3 in between the atoms 
and v i+ j +2 . In each step we would loose at most three atomic bonds by removing 
the atom but — except for the last move — also gain three new bonds by putting the 
atom back in. In the last step, however we would loose three and gain four new 
bonds, thus contradicting maximality of the number of bonds. (Since j N, we 
do not have to move the atoms . . . , Vi + j +2 themselves.) 

(ii) If <fi = I and <pj = -y for some j ^ {i — l,i + 1}, then moving the atom 
Vi to the point of £ adjacent to i>j_i, Vj and Vj + i gives three new bonds while only 
destroying two old ones, in contradiction to maximality. 

(iii) If Lfi = Lfj = I for i 7^ j, by (ii) we have (pk G {|, y , tt} for all k ^ — + 
1} U {j — 1, j + 1 } . This implies that for at most six of these k the angle (pk can be less 
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than 7r, because ^.(vi +2 —v i+ i, fj_ 2 — f <£(u,- +2 — f j+i, < ^ an d P forms 

a simple closed polygon. We therefore find some fco ^ {i — 1, i, i + 1} U {j — 1, j, j + 1} 
with tpk = ir. Now we can move the atoms Uj and to the two points of £ \ <Sjv 
adjacent to Vk with a net gain of one atomic bond. 

(iv) It is easy to see that 7^ |. An argument as in the proof of (iii) 

shows that there is a A; such that <pk = (p^+i = ir. If (pi-i = (fii+i = 7r, then 
Vi-±,Vi,Vi+i could be moved to the three points of £ \ Sn that are adjacent to 1^ 
or Vk+i- This would break seven atomic bonds but create eight new ones. Now if 
without loss of generality <fi + i is assumed to be ^r, we similarly obtain 5 — 4 = 1 
new bond by moving of Vi and v i+ i to the two points of £ \ that are adjacent to 
v k - □ 

If there is an i with cpi — |, it follows from Lemma [2.21 that deleting the atom 
Vi we are left with a graph whose boundary polygon P is in fact a convex hexagon. 
In case there is no such i, we define P to be the smallest convex hexagon whose 
segments are aligned with the bonds in £ and that contains Sn, i-e., the interior of 
P is the intersection of all half planes {x : x ■ R^i+iy^ei < c}, i = 0, . . . , 5, c G R, 

6 

containing iS^y, see Figure HI (Here Re denotes the rotation about the angle 6.) 

Let Ai, . . . , Aq G £ be the corners of P such that the segments Oj := [Ai, A i+ \] 
are numbered such that (A4+1 — Ai) ■ Ri^ei = \Ai+\ —Ai\ (mod 6). Denote the points 

of £ inside or on P by T and recall the connectedness definition from PQ: a finite 
set S C M 2 is called connected if for any two x,y G S there exist xq, . . . ,xn G S 
such that x = x, x^ = y, and the distance between successive points Xj lies 
within the interaction range of the potential, i.e. in our case \xj — Xj-i\ = 1 for all 
j- 
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Figure 4: Configuration with surrounding hexagon. 

Lemma 2.3. Let Sn be a ground state with ipi ^ | for all i. 

(i) There are points Pi, Pi G ai fl £ such that the segments {Ai, Pi) and (P/, A i+ i) 
are disjoint and 

s N nai = £n[Pi,P(\. 
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(ii) There are at most six connected components ofT\Sj^, each of them containing 
precisely one corner Ai. 

(in) Label the set of atoms in these components accordingly with Ci and denote their 
respective boundary polygons by Ci, i G {1,...,6}. Then Ci is given by the 
two segments [P/_!, Ai], [Ai, Pi] and a 'staircase' between Pi and P/_ 1; i-e., a 
polygonal path whose segments are parallel to aj_i or ai. 

(iv) #(& U . . . U C 6 ) < \A i+1 - Ai\ for all i = 1, . . . , 6. 

Proof, (i) Let Pi, resp. P[, be the point in aj HSn that is closest to Ai, resp. A i+ i. 
Then in fact Sn H a, = £ fl [Pi, P(] for otherwise there existed lattice points Vi, Vi + j 
on the polygonal path P between P, and P[ with = 4^, = . . . = V'i+j-i — 
7r,ipi + j = y> contradicting Lemma l2T27 i). 

(ii) This is a direct consequence of (i) and Theorem l2.iT ii). 

(hi) If d ^ 0, then P[_ x ^ A ^ p. Let P/_ x = Vj ,Pi = v j+r . Then ifj = 
Y' — <£j+r and for all k = j, j + 1, . . . , j + r we have (fik G n, 4^}. Let h be 
the subsequence of those indices for which cp^ ^ it. Then the angles tp^ alternate 
between ^ and 4^, i.e., tp^ = ^ -<=r- <p>k l+l = 4p for otherwise we could find 
I such that ifk t — y ~ fh+i m contradiction to Lemma I2.2( i). (Observe that 
#{/ : ifki = if} = : fk t = ^} + 1-) This implies the claim. 

(iv) If #(Ci U . . . U Cq) > \A i+ i — Ai\ + 1 = #(a« fl £), then we can successively 
move the atoms of our given ground state that lie on the segment <2j into the set 
(CiU. . .UC 6 )\(aifl£). This is possible since the number of those atoms is #((ajfl£)\ 
[C\ U . . . UC 6 )), which by assumption is not greater than #((Ci U . . . U C 6 ) \ (<2j fl £)). 
By (iii) these moves can be done in such a way that we gain three new atomic bonds 
in each step. Choosing an atom with the least neighbors among the atoms on aj, 
we will also loose only three atomic bonds in every step, except for the last one. In 
the last step, however, i.e., if there is just one atom left to move, we will only loose 
two atomic bonds, thus contradicting that Sjy was a ground state. □ 

The ground state will in general not be unique, in particular, the sets Ci in- 
troduced in the previous lemma may not consist of a single row of atoms. How- 
ever, from any ground state Sn we can always construct a 'normalized' ground 
state <Sjy by moving as many atoms as possible from P to (C\ U . . . U Cq) \ P 
as described in the proof of Lemma I2.3( iv). Moving the exterior atoms lying 
on P, we may in addition assume that in fact C 2 = ■ ■ ■ = Cq = and P = 
[A t + e 2 , A 2 ] U [A 2 , A 3 ] U . . . U [A 5 , A 6 ) U [A 6 , P 6 '] U [1%, P^ + e 2 ] U [P£ + e 2 , A 1 + e 2 ). 

We now quantify the deviation between the empirical measure of a finite-TV 
ground state and its asymptotic hexagonal shape (see (TjJ). 

This deviation can be quantified e.g. in terms of the flat norm 

||/i|| := sup |y ipdfi : Lipschitz with \<p\ < 1 and Lipschitz constant < 1^ . (7) 

Note that alternatively, instead of measuring the rate of convergence of the in 
terms the flat norm, we could start from a related absolute continuous measure (as 
used for various purposes in [1]) such as 

u N = N' 1 V XVt{x) 

' area(K(x))' 
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where V*(x) is the Voronoi cell of the point x with respect to the re-scaled triangular 
lattice N~ l / 2 C. We could then measure deviations to the limiting function /x in terms 
of the standard L 1 -distance. Due to the fact that in a ground state the rescaled 
atomic positions lie on the rate of convergence of /Xjv to /x cannot be faster 

than A r ~ 1 . On the other hand we have 

III/Jjv — HI — ||Hv — HU 1 ! < CiV -1 . 

So both approaches will yield equivalent results. 

In passing from <Sjy to Sff we have moved no more than CN^ atoms, so 

Wm-^W <CN~*, (8) 

where /x^ denotes the rescaled empirical measure of the normalized ground state 
configuration «Sjy . In fact, fluctuations of the order are to be expected due 
to surface rearrangements. In the following we will see that the optimal rate of 
convergence is, somewhat surprisingly, only of order N~i. We begin by deriving a 
corresponding upper bound. 

Theorem 2.4. Suppose Sn C £ is a ground state. There is a constant C , indepen- 
dent of N, such that /xjv, after a suitable translation, satisfies 

||a*jv - HI ^ CN~^. 

Proof. By (jHJ), without loss of generality assume that Sn = Sff- Suppose that a«, 

(n) 

i G {1, . . . , 6}, and a^, i G {1, . . . , 5} satisfy |a»/ | — |<H — 12. By construction of S N , 
ay is such that <v fl£ = DiSat if 2 < i' < 5 and [Ai + ea, A 2 ] fl£ = [Ai + e2, A 2 ] n«Sjv 
if z' = 1. We can then generate a different ground state in the following way: 

Let k G N such that |aj| + k < |dj'| — k — 1. Detach a portion W from iS^r (as 
sketched in Figure ED consisting of all the atoms in Sn which can be connected by 
at most k — 1 lattice edges of C to some lattice point on a,. Now attach U — suitably 
translated and rotated — to a^i (resp. [A\ + e2,A2] if %' = 1) restoring all atomic 
bonds. Then also the new configuration Sn is a ground state. 




Figure 5: Moving a large set of atoms keeping the number of bonds fixed. 
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Note that the union CjUC i+ i of the corresponding connected components of T\Sn 
(cf. Lemma [273]) contains at least — k— \ ai\)k lattice sites. Since |a,/| — \ at\ > 12, 
we may assume that 

-(M - KI) < fc < -(|Oi/| - |Oi|), 

so that 

#(C~ UC i+ i) > ~(\a if \ - \ ai \) 2 . 
o 

On the other hand, by Lemma l2.3( iv) and Theorem 12.1 ( hi) we know that #(C» U 
Ci+\) < C\f~N. It follows that |oj/| — |oj| < CiV* (which of course also holds true 
if \di'\ — |aj| < 12). Now using that a\ U . . . U a§ is a closed hexagon, we find that 
indeed for all i, i' G {1, . . . , 6} 

IK'I — < CN^. 
Theorem 12. If hi) now implies that 

\<k\ = )Jj + 0(N*), z = l,. ..,6. 

Rescaling by iV _ 5 concludes the proof. □ 

In order to prove by example that the bound obtained in Theorem 12.41 is sharp, 
we start by noting that the following natural candidate for an energy minimizing 
configuration is indeed a ground state. For k G N consider the regular hexagon 
spanned by B$ = B^ k \ where B\ — k{e\ — e 2 ), B 2 — fcei, £> 3 = ke 2 and B A = —Bi, 
B 5 = -B 2 , B % = -B 3 . Choose fceN such that 

3k(k + 1) + 1 < N < 3(fc + l)(fc + 2) + I. 

Now let S' N C coiw(B[ k+1 \ . . . , Bg k+1 ^) f] Che the configuration with N atoms that 
satisfies 

S' N n convfsf ] , . . . , B {k) ) = C n conv(fl? } , . . . , #f } ) 

and, in case 5^ \ conv(5? \ . . . , } ) ^ 0, #f } + ei G S' N , B {k) + e 1 -e 2 £S' N and 
the set of atoms S' N \ conv(B[ k \ . . . , Bg) is connected (see Figure [6]). 
The following Proposition is proved in |3]. 

Proposition 2.5 ([!]). a ground state. 

Note that in fact by recourse to Theorem I2.1( i) the proof of this result simply 
amounts to counting the number of atomic bonds in S' N . 

We will now show that the estimate on the rate of convergence in Theorem 12.41 
is optimal. 

Theorem 2.6. There exists a constant c > 0, independent of N, such that for 
infinitely many N G N there are two ground states S' N and Sn with N atoms such 
that 

inf {H/4 - fl N (R-+a)\\ : R G 0(2), a G R 2 } > cAT*, 
where fjf N = ± £ xe5 , 8*, fi N = ± J2 x es N 5 -t= ■ 
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Figure 6: Natural candidate for a discrete ground state. Black atoms belong to 



S' N n conv(£?j fc) , . . . , Bq K) ), white atoms belong to S' N \ conv 



>(*) 



(*) 



Proof. Let S' N be the ground state constructed above (cf. Proposition 12.51) and 
suppose iV = 3k(k + 1) + 2 for some k G N. Then 



S' N = (C n convOBi, . . . , S 6 )) U + e x } 



with = B^ 1 as above. Note that for H := h k = conv (Si, . . . , B 6 ) we have 



— ~/=XiV-V2H 



Now let m := 



and move the atoms at positions 

B 2 - se 2 + t(e 2 - ei), s,i 6 {0,...,m - 1} 

successively to the points Bi + e± + ie 2 , i — 1,2,..., m 2 in such a way that the net 
loss of atomic bonds is zero. 

Denoting the new configuration and the corresponding rescaled empirical mea- 
sure by S' N and pf N , respectively, we thus obtain that 

||AiV-/4ll < CN' 1 *. 

Now define by detaching the whole set of atoms with positions in 

{B 2 - se 2 + t(e 2 - ei) : s 6 {0, . . . , m - 1}, t 6 N} 

from S' N and attaching it — after a suitable translation and rotation — to the atoms 
on the segment [Bq,Bi] (cf. Figure EJ). 

This can be done such that the number of atomic bonds is preserved, so that we 
arrive at yet another ground state with corresponding rescaled empirical measure 
Jx N . 

If H denotes the hexagonal set 



H := conv(i?i — me 2 , B 2 — me 2 , B 3 — mei, B 4 , B 5 , B e + m[e\ — e 2 )), 



then 



Mat - ~/=Xjv-V2H 



< CN~*. 
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Figure 7: Moving a large set of atoms keeping the number of bonds fixed. 



Note that each boundary segment of H and H is of length \ ^ + 0(m), m ~ iV*. 



However, while all the segments of dH are of equal length, those of OH differ in 
length by an amount > cm > cN±. It is not hard to see that indeed 

mi{\HA(RH + a)\ : R e 0(2), a G R 2 } > civl 

(Note e.g. that diam(iJ) — diam(iJ) > y.) The assertion of the theorem now follows 
by rescaling. □ 

Note that this degeneracy of the ground state energy is in sharp contrast to the 
situation with one atom less: 

Proposition 2.7. Suppose N = 3k(k + 1) + 1 for some k. The ground state is (up 
to rigid motions) uniquely given by 

S N = Cnconv(B[ k \...,B { 6 k) ). 

So as iV grows, the situation changes rapidly between highly degenerate and 
non-degenerate energy levels. 

Proof. Suppose Sn C C is another ground state with N atoms and denote the 
corresponding set of atomic bonds by B, its boundary polygon by P. Let F and F 
be the interiors of P and P, respectively. The configurations Sn and Sn introduce a 
triangulation of F and F, respectively, into np, respectively, np, triangles of volume 
Since j^B = j^B and, by Theorem 12.11 F and F are simply connected, we have 
by Euler's formula uf = np and hence \F\ = \F\. Again by Theorem 12.11 we have 

< H 1 (P) = U 1 (P). 

We now use the following three facts. First, by inspection the above perimeters 
can be expressed in terms of the following surface energy functional derived in [1] 
as a Gamma-limit of the Heitmann-Radin energy: 



2V}(P) = J T(v), 2V}{P) = J T(u) 



9F dF 
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where v denotes the outward unit normal to the domain of integration and T is the 
7r/3-periodic function defined by 



T(u) = 2(z/ 2 - vJy/3) for v 



( 



COS (f 



sin if 



) 



v e [o,tt/3]. 



(9) 



Second, such surface energy functionals have a unique minimizer given by the Wulff 
construction, as proved by Taylor (in the language of geometric measure theory) 
and Fonseca-Muller (in the present language of boundary integrals): 

Theorem 2.8. (HI 12] A functional of form 



with T : S*™ -1 — > [0, oo) continuous and bounded away from zero, is minimized over 
sets E C W 1 of finite perimeter and volume c > if and only if E agrees, up to 
translation and up to a set of measure zero, with \ c Wr, where Wr is the Wulff set 



and X c > is the unique normalization constant such that XWr has volume c. 

Finally, for the energy (Q, an elementary calculation shows that the Wulff set 
is given by the intersection of the six half-spaces x ■ v < T(u) for the minimizing 
normals ^(-f^), j = 0, . . . , 5. 

Since F = \\f\Wt, this establishes that both F and F are minimizers of / 
subject to the same fixed volume, so that by Theorem 12.81 they coincide and therefore 
Sn = Sn (up to translation). □ 
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